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We report on the investigation of an advanced circular plasmonic nanoantenna under ultrafast 
excitation using nonlinear photoemission electron microscopy (PEEM) under near-normal incidence. 
The circular nanoantenna is enhanced in its performance by a supporting grating and milled out 
from a gold film. The considered antenna shows a sophisticated physical resonance behavior that 
is ideal to demonstrate the possibilities of PEEM for the experimental investigations of plasmonic 
effects on the nanoscale. Field profiles of the antenna resonance for both possible linear polarizations 
of the incident field are measured with high spatial resolution. In addition, outward propagating 
Hankel plasmons, which are also excited by the structure, are measured and analyzed. We compare 
our findings to measurements of an isolated plasmonic nanodisc resonator and scanning near-field 
optical microscopy (SNOM) measurements of both structures. All results are in very good agreement 
with numerical simulations as well as analytial models that are also discussed in our paper. 


I. INTRODUCTION 

Optical nanoantennas provide a possibility to con¬ 
vert light between the far-field and nanoscale dimensions 
and constitute building-blocks of utmost importance in 
nanooptics [1-4]. They are typically created from sub¬ 
wavelength plasmonic resonators. A variety of different 
design approaches has emerged, from nanowire [5] or 
bow-tie antennas [6] to sophisticated geometries where 
the emission characteristic and resonance strength can 
be tailored by design [7]. The short response time of 
the involved plasmonic excitations allows to enter the 
realm of ultrafast light control at the nanoscale with these 
structures, opening entirely new research perspectives for 
physics, chemistry, biology, and material science [8, 9]. 

The success of such approaches critically depends on 
the possibility to investigate the structures experimen¬ 
tally on the nanoscale. Scanning near-field optical micro¬ 
copy (SNOM) has become an established and accepted 
method for this purpose [10, 11]. However, the presence 
of a physical probe, which disturbs the optical function¬ 
ality of the structure and limits the spatial resolution, as 
well as the scanning nature of the method prompted the 
development of alternative experimental techniques [12]. 

Photoemission electron microscopy (PEEM) has be¬ 
come a major tool for this purpose in recent years when 
combined with ultrafast light sources [8, 9, 13-20]. The 
interference of surface plasmons and the incident light 
pulse allows the direct visualization of the plasmon field 
in the nonlinear regime, i.e. when multiple photons 
are required to excite one photoelectron. Since most 
PEEM experiments are done under a grazing angle of in¬ 
cidence, wide beating patterns and Moire-like images are 
observed when imaging the straight propagation of sur¬ 


face plasmons excited by slits, trenches, or other straight 
edges [17, 20-24]. The interference from structures with 
circular symmetry like nanoholes, however, show wide 
parabolic ring-like patterns [25, 26]. This limitation on 
the achievable image quality is mitigated when normal 
or near-normal incidence is used [27]. 

The aforementioned publications concentrated on the 
investigation of simple propagating surface plasmon ge¬ 
ometries or hot-spots generated by nanoparticles. The 
study of more complex systems with a certain nanoopti- 
cal functionality using nonlinear PEEM has so far been 
limited to the investigation of bow-tie nanoantennas [28] , 
plasmonic focusing devices [29], nanoparticle-on-plane 
geometries [30], or cross resonant optical antennas [31]. 

In this work, we wish to investigate the behavior of 
a circular optical nanoantenna under ultrafast excitation 
using nonlinear PEEM at near-normal incidence. The 
resonance of the antenna was enhanced by a properly 
designed surrounding grating [7]. We find a pronounced 
enhancement of the photoelectron yield both for TE and 
TM polarization of the incident light. Our results are 
compared to PEEM measurements on an isolated nan¬ 
odisc antenna as well as SNOM measurements on both 
structures. PEEM is found to be able to resolve the 
field-profile of the antenna modes with much better spa¬ 
tial resolution than SNOM. The ring-like structure also 
supports the excitation of outward propagating Hankel 
plasmons on the surrounding gold film [32], which we 
also observe for both polarizations. The measurements 
are accompanied by rigorous numerical simulations as 
well as analytical modeling of the interference of Hankel 
plasmons and the illumination under ultrafast excitation 
conditions, which explain the experimental findings with 
very good agreement. 
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Figure 1. Scheme of the used sample and experimental con¬ 
figuration. A central nanodisc is surrounded by a circular 
grating with period p. The grating helps enhancing the reso¬ 
nance in the central disc. The entire structure is surrounded 
by a gold him and the configuration is probed by an ultrashort 
laser pulse. 


II. CIRCULAR OPTICAL NANOANTENNA 
WITH ENHANCED RESONANCE 

We start by introducing the physical working principle 
of the considered structure. A plasmonic nanoantenna 
generally consists of a small metallic nanostructure that 
can sustain a plasmonic mode which is able to oscillate in 
resonance. The resonance frequencies and other proper¬ 
ties such as the held prohle or emissivity of the nanoan¬ 
tenna are determined by the size, shape and material of 
the structure and can be tuned to ht a particular intended 
functionality [2]. 


(1 2 ) 

second kind and order l are denoted by ’ , where l is 
an integer. With a mathematical approximation for the 
Hankel function [35], an expression for the total electric 
held of the HSPP is found 

E hspp (r,w) = [fc hspp r] -5 e* fchsp p r fe z + i-^—e r ^\ , 

\ ^hspp J 

(3) 

where k = [^ spp — . Vectors represent unit 

vectors in i-direction. We have neglected an additional 
term 

*tr& + a ") <4> 

since it is quickly damped away with growing distance 
from the center (hrst several oscillations). We may thus 
call (3) a ’’far-held” approximation of the Hankel plas- 
mon. 

The localized resonance of the nanodisc can be under¬ 
stood by regarding the structure as a modal Fabry-Perot 
resonator [36]. The outward-propagating HSPPs are par¬ 
tially reflected back into inward-propagating modes at 
the circumference. This modal reflectivity can be tailored 
by proper design of the mode mismatch at the circum¬ 
ference [37]. The resonance condition is a constructive 
interference and reads 

k hspP R + $ = ^ • 7r, (5) 

where R is the radius of the disc and <F is the phase 
acquired upon reflection at the circumference. The radial 
mode number is denoted as v and needs to be an integer. 
We use v = 1 in our case. 


A. Central gold nanodisc 


B. Surrounding ring grating for enhanced 
resonance 


Our basic design is shown in Fig. 1. The central res¬ 
onant element is a circular gold nanodisc. For such a 
system with circular symmetry, in- and outward prop¬ 
agating Hankel surface plasmon polaritons (HSPPs) are 
the fundamental solutions for the electromagnetic surface 
modes. At a hxed frequency, the electric field component 
perpendicular to the surface is described as [33] 

E z { r, w) oc A 1>2) [ fc hsppHr] e ll,fi , (1) 

where r = [x 2 + y 2 ] i is the radial coordinate and cp de¬ 
notes the angular coordinate in the x-7/-plane.The prop¬ 
agation constant of the plasmon is denoted as & hspp and 
reads as [33, 34] 


& hsPP M = 


<bnM 2 

^m(^) + 1_ 


( 2 ) 


where e m is the permitivity of gold. Vacuum was as¬ 
sumed as cladding. The u -dependence of k hspp is not ex¬ 
plicitly noted hereafter. The Hankel functions of first or 


The resonance strength of a nanoantenna can be in¬ 
creased significantly when the central element is enclosed 
by a periodic structure [7, 38]. The configuration acts as 
a Bragg reflector if concentric rings are used in our case 
to enclose the central nanodisc. One can tune the peri¬ 
odicity such that the condition for resonant coupling of 
outward- into inward-propagating HSPPs is met. This 
Bragg condition for HSPPs reads as 

KsppP = m-n, (6) 

where m is an integer denoting the band index and p 
is the period of the grating. The grating increases the 
resonance strength in the central disc by enhancing the 
partial reflection of outward-propagating HSPPs at its 
circumference. This increase in modal reflectivity en¬ 
hances the quality factor of the resonator and leads to a 
better performance of the nanoantenna. 

A second possibility to employ the grating considers 
the excitation of the antenna by incident light. The bare 
disc can be excited from free-space illumination in the 
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quasi-static regime thanks to its subwavelength size [1]. 
Such an approach does not require phase matching to the 
Hankel modes and can be used to excite the disc at any 
wavelength within the bandwidth of its localized plas- 
monic resonance. However, it is also possible to employ 
the grating to couple incident light resonantly into Han¬ 
kel plasmons [7]. The prerequisite for this is known as 
the Bloch condition and reads as 

27r 

k 0 smO + n -= & hspp , (7) 

P 

where n is an integer. The angle of incidence is denoted 
by 0. 

Our design approach originates from the simple obser¬ 
vation that (6) and (7) can be fulfilled simultaneously for 
certain combinations of 0 , m, and n if 

/ 2 n\ 

k hspp ( 1 - — J - k 0 sm6 = 0. (8) 

Since we want to excite the structure from near-normal 
incidence (0 ~ 0°), we use a grating corresponding to 
m = 2 and n = 1. The period in our structure was 
thus chosen to be equal to the HSPP wavelength p = 

27r//c hspp ^hspp • 

In this way, the circular grating works as a three-mode- 
coupler. The occurrence of a Bragg resonance facilitates 
the interaction of the inward- and outward-propagating 
Hankel mode which manifests itself by the formation of a 
bandgap. The HSPPs become evanescent in the grating 
region and a coupling between the two otherwise orthogo¬ 
nal modes is possible within the penetration depth of the 
modes into the grating. While the usually used first order 
Bragg resonance is located at the band edges, the second 
order Bragg resonance appears at the T-point, i.e. for 
a vanishing in-plane quasi-wavevector. This facilitates 
the interaction of both the in- and outward HSPP mode 
with near-normal incident light. The coupling to inward 
propagating HSPPs will further enhance the resonance 
of the nanoantenna whereas the excited outward propa¬ 
gating HSPP will be visible on the gold film surrounding 
the structure. 


III. PEEM MEASUREMENTS 

A. Sample preparation and experimental setup 

We investigated the nanoantenna structure by means 
of multiphoton photoemission electron microscopy 
(PEEM). Samples for the ring-enhanced nanoantenna 
structure and also the bare nanodisc antenna were fab¬ 
ricated by evaporating a gold film on a fused silica sub¬ 
strate. The film thickness exceeded 200 nm to exclude 
effects from the substrate side for this sample. The struc¬ 
tures were milled into the gold film using Focused Ion 
Beam Milling (FIB) with gallium ion irradiation. The 
width of the milled grooves was 80 nm, the central disc 
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Figure 2. Sketch of the PEEM setup. A polarized 120 fs 
pulse is focused on the sample S with a lens L (focal length 
250 mm). The polarization is controlled by a half-wave plate 
(HWP). The beam was stabilized using the beam stabiliza¬ 
tion system Aligna (TEM Messtechnik Germany), allowing 
for stability during long-term measurements. A part of the 
beamline was passed to a beam position detector (BPD) be¬ 
hind a beam splitter (BS) for this purpose and the feedback 
was used to control the motorized mirrors (MM). The chirped 
mirrors (CM) compressed the pulse to the final duration of 
30 fs. 


had a diameter of 180 nm and the period of the grating 
was 780 nm. 

The experimental setup is shown in Fig. 2. As 
laser source we used a mode-locked Ti:Sa oscillator 
(Coherent Micra), which provides an average output 
power of 100 mW at the sample position at a repetition 
rate of 80 MHz, resulting in a pulse energy of approxi¬ 
mately 1 nJ. The duration of the ultrashort pulses was 
30 fs, the central wavelength was 800 nm and the spectral 
bandwidth 60 nm. 

The frequency of the laser oscillator corresponds to a 
photon energy of 1.55 eV. The work function for gold 
evaporated in our facilities is in the order of 4.6 eV. This 
means that photoemission of electrons is just possible 
in a nonlinear 3-photon process which requires enough 
intensity in the optical field in order to be observed. 

We used a PEEM from Focus GmbH (Germany) to 
image the emitted photoelectrons. The illumination an¬ 
gle of incidence was 4° with respect to the sample nor¬ 
mal. A wide Gaussian beam was imaged onto the sample 
using a / = 250 mm lens, resembling a quasi-plane-wave 
illumination on the sample. Photoexcited electrons emit¬ 
ted from the sample are accelerated by a static extractor 
voltage up to 16 kV and imaged by an electrostatic lens 
system which provides up to 20 nm lateral resolution de¬ 
pending on the field of view. The amount of emitted 
electrons can be expected to be small enough in order to 
not influence each other. The total exposure time was 
about 15 min for the ring-enhanced structure. This was 
optimized to give the best signal-to-noise ratio in the ex¬ 
periments. 
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Figure 3. Logarithmic plot of the PEEM yield for the ring- 
enhanced sample in (a) TM and (b) TE polarization. Both 
polarizations show noticeable yield but only TM polarization 
excites the central disc. The dashed lines denote the areas 
from which the data for the detailed analysis of outward- 
propagating Hankel plasmons visible in the outer parts was 
extracted (Sec. HID). 


B. Results and interpretation of the PEEM 
measurements 

The experimental findings for the ring-enhanced struc¬ 
ture in TE and TM polarization are shown in Fig. 3. 
We observe first of all a dark area around the sample 
structure, which is caused by the SEM inspection dur¬ 
ing fabrication. The scanning electron beam provides 
the energy to permanently adsorb chemicals (most likely 
carbon) which are present. Since PEEM is extremely 
sensitive to surface conditions, even a monolayer of ad¬ 
sorbed material can alter the effective work function. In 
our case, it is slightly increased in the scanned area which 
explains the darker region in the center. The outer region 
represents parts which were not scanned by the electron 
beam. 

Although the effect seems disturbing it actually in¬ 
creases the signal-to-noise ratio of our measurements and 
provides the possibility to accurately overlay PEEM im¬ 
ages with the SEM image. This allows us to create the 


Figure 4. PEEM images overlayed with the SEM image for 
(a) TM and (b) TE polarization. A clear qualitative dif¬ 
ference in the field localization and excitation of the central 
nanodisc is observed as well as details of the localization of 
the photoemission. 


overlay plots shown in Fig. 4, where the origin of the 
PEEM yield can be studied in great detail with high spa¬ 
tial resolution. We estimate the positioning accuracy of 
the overlay to be about 60 nm. 

Photoemission is known to be very sensitive to the 
presence of an electric field component perpendicular to 
the surface [13]. In the case of TM illumination, we thus 
find a strong excitation of the central disc. Both polar¬ 
izations were ensured to carry the same power and were 
measured using the same exposure time. Apart from the 
strong resonance in the center, we also find a weaker lobe¬ 
like enhancement in the horizontal direction, which was 
the plane of incidence. This is consistent with Hankel 
plasmon theory, which predicts that the excitation in¬ 
tensity drops as cos 2 ((/?) for TM modes, where (p is the 
angle to the horizontal direction. Consequently, no ex¬ 
citation of HSPPs is observed in the vertical direction. 
It should also be noticed that the photoelectron yield is 
slightly increased at the positions of the rings. Although 
there cannot be electron emission from the gaps, sharp 
edges lead to a higher probability for electron emission 
due to confinement effects. This manifests itself in the 
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increased yield at the ring positions since the gaps are 
too narrow to be accurately resolved. 

We were able to measure a weaker but noticeable yield 
also for TE polarization. The reason for this is, that the 
excited HSPPs themselves provide for the necessary per¬ 
pendicular electric field component. This argumentation 
is supported by the fact that we measured a much lower 
yield for the bare nanodisc sample in TE polarization 
as well, where the excitation of HSPPs will be shown to 
be only very weak in Section HIE. In the ring-enhanced 
case, a clear excitation along the vertical image direction 
is visible. For this polarization, theory predicts a sin 2 ((^) 
dependence of the intensity, which is in agreement with 
our experiments. However, we find that the picture has 
clear qualitative differences to the TM case. The central 
disc was not excited. Instead, a larger dipolar Hankel 
mode in vertical direction forms on the area between the 
disc and the first ring. We again observe a lobe-like spa¬ 
tial structure in polarization direction and an increased 
yield from the ring positions, similar to the TM case. 


action Hamiltonian p • E, where p is an atomic dipolar 
moment and thus leaves us with the overall dependence 
on the electric field to the sixth power in the 3-photon 
case. This approach differs from using the third power 
of the local intensity, since intensity is defined as a time- 
average of the Poynting vector [41] 


I = 


i R e 


E(r,cj) x H*(r ,c 


(9) 


although the principle overall power dependence Y oc I 3 
stays intact and differences occur just in the weights of 
certain polarization components contributing locally to 
the PEEM yield image. 

The total value of the instantaneous electric field in¬ 
volves polarization contributions perpendicular and par¬ 
allel to the sample surface 


|E(r, t)\ 


Ej_(r ,t) +£f(r,f) 


( 10 ) 


We therefore arrive at our used model formula for the 
C. Comparison to FDTD data 3-photon PEEM yield 


We compared the experimental PEEM images to rig¬ 
orous simulations using the finite-difference time-domain 
method (FDTD). Full 3D simulations were employed un¬ 
der 4° incidence using the freely available FDTD software 
meep [39]. The simulation time was chosen long enough 
for the field to reach the steady-state. The results are 
shown in Fig. 5. 

The left panel shows the intensity calculated close to 
the sample surface for both polarizations. It becomes ev¬ 
ident that the qualitative difference we observed between 
the two polarization directions is inherent to the inves¬ 
tigated structure. The TM simulation confirms that the 
central disc is the location of the highest intensity in this 
case whereas it is not excited in the TE case. The sym¬ 
metry of the system is broken by the 4° incidence angle 
while the additional rings increase the angular sensitivity 
of the antenna. 

Although the FDTD intensity results confirm in princi¬ 
ple our PEEM measurement data, the images look quan¬ 
titatively different. The background level within the ring 
area is much higher in the measurement and no pro¬ 
nounced two-lobe structure is seen for the TE case. To 
explain this discrepancy, one must not forget that the 
photoelectron yield is not a direct measure for the elec¬ 
tromagnetic fields. An accurate modeling is needed to 
connect the two quantities. 

We therefore employ a model which connects the pho¬ 
toelectron yield to the temporal integral of the sixth 
power of the instantaneous electric field. This is facili¬ 
tated by the fact that the probability for three-photon 
photoemission should be proportional to the absolute 
value squared of the matrix element for three photon ab¬ 
sorption, following a nonlinear generalization of Fermi’s 
golden rule [40]. This matrix element involves the inter- 


Y( r)« J |E(r, t)| 6 di = J [^(r, t) + £fi (r, t)]* dt. 

( 11 ) 

For the Hankel plasmon, both polarization contribu¬ 
tions are present in form of a perpendicular and radial 
component and approximately 7r/2 out of phase. The 
illuminating plane wave either contains only a parallel 
component (TE case), or both components in phase (TM 
case) which are superimposed during the duration of the 
pulse. Since we use just a single ultrashort 30 fs pulse 
to perform the experiment, HSPPs are excited from the 
early parts of the pulse and interfere with later parts 
in a spatio-temporal coherent manner. This creates the 
interference patterns typical for multiphoton PEEM ex¬ 
periments. 

To take this into account, we modeled the PEEM yield 
by coherently superimposing the nanoantenna near-field 
extracted from the FDTD data with the corresponding il¬ 
luminating field components, allowing the amplitude and 
total phase of the latter as a free fit parameter. 

The last panel of Fig. 5 shows the result of the model¬ 
ing. The same fit parameters were used for both polar¬ 
izations. A remarkable agreement with the measurement 
data is achieved. The TE case shows the correct two-lobe 
structure now. Additionally, background levels and even 
finer details are accurately reproduced qualitatively and 
quantitatively. 

A peculiar difference is that the model overestimates 
the yield from the positions of the rings. The reason 
for this is the plasmonic field enhancement within small 
gaps, which is present in the FDTD data. Since the gaps 
are not filled with metal, no enhancement of the photo¬ 
electron yield can occur at these positions apart from the 
sharp edges as discussed earlier. 
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I 3 (from FDTD) PEEM measurement yield model (from FDTD) 


Figure 5. Comparison between the third power of the intensity from FDTD data, the actual PEEM measurement and the yield 
computation from the model described in the text. The colorbars apply to all plots in the same row. The yield model is in 
very good agreement with the measurement, whereas the intensity data does not give the correct background levels. 


D. Analysis of propagating Hankel plasmons 

In addition to the antenna resonance, we observe 
outward-propagating HSPPs on the surrounding gold 
film in both polarization cases. They can be analyzed 
in detail since a full analytical theory is readily available 
[33]. We confirmed that the measured HSPPs obeyed the 
predicted angular dependence of the intensity as cos 2 ((^) 
and sin 2 ((/?), respectively (not shown). 

To further investigate their properties, we extracted 
a cut through the horizontal (TM) or vertical direction 
(TE), shown as dashed areas in Fig. 3. We plot the results 
in Fig. 6. In the TE case, the HSPPs propagate perpen¬ 
dicular to the plane of the incident pulse. The measured 
oscillation period therefore matches exactly the plasmon 
wavelength A hspp . For TM polarization, the HSPP pulse 
propagates parallel or antiparallel to the excitation pulse 
direction, respectively. The wavenumber of the oscilla¬ 
tion is thus altered by ±kosin6 in this case. However, 
this effect is only minimal, since sin(4°) « 0.07. 

The attenuation of the oscillations is governed by two 
terms. The first contribution comes from the physics 
of the Hankel plasmon itself. As radial solutions, they 
drop with a factor [k hspp r ]~2 quicker than plane surface 
plasmon polaritons. The second factor originates from 
the limited illumination time during the ultrashort pulse. 
If we assume a Gaussian envelope, the temporal width 
r of the pulse translates into a spatial envelope for the 
excitation which leaves us with the total dependence of 



Figure 6. Detailed analysis of outward propagating Hankel 
plasmons from the regions marked in Fig. 3 and comparison to 
the analytical model. The yield envelope contains information 
about the temporal structure of the pulse as well as the spatial 
decay of Hankel plasmons. 
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the envelope 


^envO) OC 


R e(fc h 


exp 



■ const. 


( 12 ) 

where n g is the group index of the plasmon pulse and r 
is the distance to the origin. The vacuum speed of light 
is denoted by c. From the dispersion relation of HSPPs 
at 800 nm free-space wavelength, we calculate n g = 1.08. 
The term re is approximately 9 pm and represents a char¬ 
acteristic length scale for the experiments. The analyt¬ 
ical model fits the experimental data excellently. This 
also opens the possibility to check or even measure the 
pulse length of the excitation from the PEEM images. 

We emphasize that equation (12) does not contain 
a contribution which is inherent to the intrinsic ohmic 
metal loss, i.e. a term proportional to exp[—Im(fc hspp )r]. 
This term can be neglected since 


Im(fe hS pp) ^ 2 q—3 
Re(&hs PP ) ~ 


(13) 


and does not contribute noticeably over the distances rel¬ 
evant in our experiment. 

We use the yield model described above to fit the ex¬ 
perimental data. The total electric field is a superposi¬ 
tion of the HSPP and the illuminating plane wave with 
excitation amplitudes A and B 


E(r, t) = A- E hspp (r, t) + B • E pw (r, t), (14) 


where both terms contain a moving Gaussian pulse en¬ 
velope 


E pw (r,t) = E pw (r,t) • exp 
E h8 pp(r,i) = E happ (r,t) • exp 


( t — sin Ox/c) 2 


(t - n g r/c ) 2 


(15) 

(16) 


in addition to their steady-state solutions 
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Figure 7. Logarithmic plot of the PEEM yield for the bare 
nanodisc sample in TM polarization. An enhancement of pho¬ 
toemission due to the plasmonic resonance is observed at the 
position of the disc which is weaker than in the ring-enhanced 
case. 



the background levels far from the center were not ideally 
equal at all positions. We believe that a slight spatial in¬ 
homogeneity of the illumination at the sample position 
is the reason for this since it is a wide Gaussian beam 
in reality and not an ideal plane wave. Some deviations 
from the measurement also occur for smaller distances 
from the center which happens since the wave originates 
from different rings in the experiment and not just from 
an ideal point source. The approximation introduced for 
HSPPs in (3) therefore does not fully hold close to the 
outermost ring. Nevertheless, we were able to show that 
the PEEM images contain rich spatio-temporal informa¬ 
tion which allows a detailed analysis and understanding 
of the underlying physics of the structure. 


E. Comparison with measurements on the bare 
nanodisc 


Epw(M) = - 

E hspp (r,t) = ^ 


e ik 0 sin Ox-iut + c c 

Eh S pp (r, w) e~ luJt + c.c. 


(17) 

(18) 


E(, denotes the polarization vector in the TE or TM case 


E 0 = 


[0,1,0] (TE) 

[cos#, 0, sin#] (TM) 


(19) 


and E hspp (r,cj) the analytical HSPP solution given in 
eq. (3). 

The yield model (11) is used with the above field ex¬ 
pressions to fit the excitation amplitudes to the experi¬ 
mental data in Fig. 6. The results are in very good agree¬ 
ment with the experimental findings. Slightly different 
fitting values had to be used in the different cases since 


The results for the ring-enhanced structure are com¬ 
pared to the measurements on the bare nanodisc in order 
to quantitatively verify the enhancement mechanism of 
the antenna. The total exposure time had to be increased 
to ~ 2 h in order to get a significant signal. Fig. 7 shows 
the experimental result for the bare plasmonic nanodisc. 
The polarization of the incident light was TM in this 
case, the electrical field vector being in the plane of in¬ 
cidence. The yield of photoelectrons at the position of 
the nanodisc was found to be much stronger than on the 
surrounding planar gold film, but much weaker than in 
the ring-enhanced case. This is caused by the localized 
quasi-static plasmonic resonance of the nanodisc, in ac¬ 
cordance with our theory. The incident plane was not 
exactly horizontal in Fig. 7 which can be seen from the 
weak excitation of HSPPs excited by the single groove 
and propagating away from the disc. However, the exci¬ 
tation strength is so low that we could not observe notice- 





















able yield for TE polarization using the same exposure 
time (not shown). In this case, the HSPP is the only 
source of a normal electric field component. 

If we compare the maximum values for E 1 / 3 in the TM 
case and take into account the different exposure times, 
we find an about 6-fold enhancement of the near-field 
intensity at the nanodisc by adding the ring grating. 


IV. COMPARISON WITH SNOM 
MEASUREMENTS 

The PEEM results are compared to previous measure¬ 
ments [7] on the same structure using scanning near¬ 
field optical microscopy (SNOM) as established and ac¬ 
cepted experimental method for nanophotonics. We had 
to adapt the initially described sample design slightly 
in order to account for the experimental circumstances. 
The experiments were performed with a diode laser of 
785 nm wavelength which illuminated the structure from 
the substrate side. A SNOM tip with a circular aperture 
was used to scan the sample from the top. The thickness 
of the gold film had to be decreased to 50 nm in order 
to allow sufficient coupling of light to the upper side. 
This results in a certain interaction between plasmons at 
the substrate and cladding interfaces. It was shown in 
a previous work, that this interaction leads to a slightly 
altered value for & hspp , but does not influence the basic 
functionality of the structure much, neither qualitatively 
nor quantitatively if we regard the field profiles and en¬ 
hancement factor, which are our values of interest here 

[7]. 

A wide and polarized Gaussian beam was imaged at 
the sample position so that the excitation resembles a 
plane wave. The aperture tip was gold coated and 
light was collected into a single mode optical fiber which 
guided the signal to a photodetector. 

Fig. 8 shows a summary of the SNOM results. 

We started by investigating the single disc first. A clear 
excitation of a localized plasmonic resonance at the disc 
position was observed, accompanied by a small amount 
of background light from direct transmission of the illu¬ 
mination. The structure surrounded by the rings showed 
a much more complex behavior. We found pronounced 
side lobes which we attribute to the excitation of Hankel 
plasmons by the grating. They appear only in horizontal 
direction in the images which is consistent with the inci¬ 
dent polarization direction and the PEEM findings. We 
also observe an enhancement of the field strength at the 
central disc. The near-field intensity was again approx¬ 
imately 6 times higher for the ring-enhanced structure 
than for the bare nanodisc. This shows that the PEEM 
measurements lead to quantitatively comparable results. 

The SNOM data was also compared to rigorous FDTD 
simulations. We had to account for certain experimen¬ 
tal circumstances by processing the raw simulation data. 
The tip is sensitive to electric field components perpen¬ 
dicular to the tip axis and approaches the sample under 


an angle of 30° [42]. Furthermore, the apex size of ap¬ 
proximately 180 nm limits the achievable resolution. This 
means that all images get slightly blurred as if they were 
taken further away from the sample. We took that into 
account by projecting the FDTD data on a plane inclined 
by 30° with respect to the surface normal and used an 
extraction distance of 100 nm above the structure. The 
results are shown in Fig. 8 and show a good agreement 
with the experimental findings. This demonstrates the 
superior possibilities of PEEM over SNOM with respect 
to resolution. In both cases, accurate modeling of the 
underlying physics is necessary to interpret the result¬ 
ing images and compare to, e.g., rigorous simulations or 
analytical data. 


V. CONCLUSION 

We investigated a circular optical nanoantenna un¬ 
der ultrafast femtosecond laser irradiation be means of 
multiphoton photoemission electron microscopy. Posi¬ 
tions at which the electromagnetic field got enhanced 
by the constructive interference of Hankel plasmons were 
found to lead to a much higher yield in photoelectrons. 
The method thus provides information about the spatial 
structure of the electromagnetic field with high accuracy 
on ultrashort timescales. 

For our investigations, we used an antenna design 
composed of a plasmonic nanodisc, surrounded by rings 
which function as a grating. It acts as a Hankel plas- 
mon coupler on the gold film used for the experimen¬ 
tal realization of the device. Inward propagating Hankel 
plasmons are resonantly excited by the illumination and 
drive a localized resonance at the central disc of the struc¬ 
ture, where the grating confines the energy by a Bragg 
resonance. This enhances the near-field intensity by a 
factor of ~ 6 compared to the localized plasmonic res¬ 
onance of the disc alone. The grating also excites out¬ 
ward propagating Hankel plasmons which were observed 
in the PEEM and shown to contain information about 
the propagating plasmons themself, as well as the ultra- 
short excitation of the structure by a femtosecond pulse. 

We compared our findings to measurements of the 
same antenna design using scanning near-field optical 
microscopy (SNOM). While this method has the advan¬ 
tage of directly mapping the electromagnetic near-field, 
it was shown that the physical presence of a probe has 
disadvantages regarding the experimental configuration 
and achievable resolution. Moreover, the measurement 
of ultrashort electromagnetic phenomena is very difficult 
due to the scanning nature of the method. The achieved 
PEEM images had a much better spatial resolution than 
the SNOM images and do not require any disturbing 
probe or sample scanning which leaves much freedom to 
the possible excitation scenarios. However, it was shown 
that accurate physical modeling is necessary to relate the 
acquired PEEM images to the electromagnetic field since 
the process is of nonlinear nature. 
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Figure 8. Comparison of SNOM measurement results and FDTD modeling for the bare nanodisc sample and the ring-enhanced 
structure. Measurement and model verify the enhancement of the near-field intensity. The scalebar is 2 pm long in each image. 
The spatial resolution is limited to 180 nm, which was the apex size of the aperture SNOM tip. 



Since the photon energy of the used 1.55 eV radiation 
is much lower than the workfunction of gold, a 3-photon 
process was necessary to excite photoelectrons. The pho¬ 
toelectron yield was modeled as the temporal integral of 
the instantaneous yield, which was shown to depend on 
the total electric field to the sixth power in our case. This 
field was composed of an electromagnetic interference be¬ 
tween both the plasmon fields and the illumination itself. 
Rigorous electromagnetic simulations were performed us¬ 
ing the finite-difference time-domain method (FDTD). 
An analytical theory was employed as well to describe the 
outward propagating Hankel plasmons. Both were found 
to be in very good agreement with the experimental im¬ 
ages when we apply our yield model to the theoretical 
and numerical data. 

Nevertheless, more work has to be done in the future 


to clarify the exact role of the local polarization as a func¬ 
tion of time for the acquisition of the multiphoton PEEM 
images. This includes the implementation of a pump- 
probe configuration in our setup to probe the complete 
spatiotemporal dynamics of the structure. 
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